Comparisons of host mitochondrial, nuclear and endosymbiont bacterial genes reveal cryptic fig wasp species and the effects of Wolbachia on host mtDNA evolution and diversity

Background Figs and fig-pollinating wasp species usually display a highly specific one-to-one association. However, more and more studies have revealed that the "one-to-one" rule has been broken. Co-pollinators have been reported, but we do not yet know how they evolve. They may evolve from insect speciation induced or facilitated by Wolbachia which can manipulate host reproduction and induce reproductive isolation. In addition, Wolbachia can affect host mitochondrial DNA evolution, because of the linkage between Wolbachia and associated mitochondrial haplotypes, and thus confound host phylogeny based on mtDNA. Previous research has shown that fig wasps have the highest incidence of Wolbachia infection in all insect taxa, and Wolbachia may have great influence on fig wasp biology. Therefore, we look forward to understanding the influence of Wolbachia on mitochondrial DNA evolution and speciation in fig wasps. Results We surveyed 76 pollinator wasp specimens from nine Ficus microcarpa trees each growing at a different location in Hainan and Fujian Provinces, China. We found that all wasps were morphologically identified as Eupristina verticillata, but diverged into three clades with 4.22-5.28% mtDNA divergence and 2.29-20.72% nuclear gene divergence. We also found very strong concordance between E. verticillata clades and Wolbachia infection status, and the predicted effects of Wolbachia on both mtDNA diversity and evolution by decreasing mitochondrial haplotypes. Conclusions Our study reveals that the pollinating wasp E. verticillata on F. microcarpa has diverged into three cryptic species, and Wolbachia may have a role in this divergence. The results also indicate that Wolbachia strains infecting E. verticillata have likely resulted in selective sweeps on host mitochondrial DNA

Wolbachia-Induced Unidirectional Cytoplasmic Incompatibility and Speciation: Mainland-Island Model

Bacteria of the genus Wolbachia are among the most common endosymbionts in the world. In many insect species these bacteria induce a sperm-egg incompatibility between the gametes of infected males and uninfected females, commonly called unidirectional cytoplasmic incompatibility (CI). It is generally believed that unidirectional CI cannot promote speciation in hosts because infection differences between populations will be unstable and subsequent gene flow will eliminate genetic differences between diverging populations. In the present study we investigate this question theoretically in a mainland-island model with migration from mainland to island. Our analysis shows that (a) the infection polymorphism is stable below a critical migration rate, (b) an (initially) uninfected “island” can better maintain divergence at a selected locus (e.g. can adapt locally) in the presence of CI, and (c) unidirectional CI selects for premating isolation in (initially) uninfected island populations if they receive migration from a Wolbachia-infected mainland. Interestingly, premating isolation is most likely to evolve if levels of incompatibility are intermediate and if either the infection causes fecundity reductions or Wolbachia transmission is incomplete. This is because under these circumstances an infection pattern with an infected mainland and a mostly uninfected island can persist in the face of comparably high migration. We present analytical results for all three findings: (a) a lower estimation of the critical migration rate in the presence of local adaptation, (b) an analytical approximation for the gene flow reduction caused by unidirectional CI, and (c) a heuristic formula describing the invasion success of mutants at a mate preference locus. These findings generally suggest that Wolbachia-induced unidirectional CI can be a factor in divergence and speciation of hosts

CrossHybDetector: detection of cross-hybridization events in DNA microarray experiments

Background\ud DNA microarrays contain thousands of different probe sequences represented on their surface. These are designed in such a way that potential cross-hybridization reactions with non-target sequences are minimized. However, given the large number of probes, the occurrence of cross hybridization events cannot be excluded. This problem can dramatically affect the data quality and cause false positive/false negative results.\ud \ud Results\ud CrossHybDetector is a software package aimed at the identification of cross-hybridization events occurred during individual array hybridization, by using the probe sequences and the array intensity values. As output, the software provides the user with a list of array spots potentially 'corrupted' and their associated p-values calculated by Monte Carlo simulations. Graphical plots are also generated, which provide a visual and global overview of the quality of the microarray experiment with respect to cross-hybridization issues.\ud \ud Conclusion\ud CrossHybDetector is implemented as a package for the statistical computing environment R and is freely available under the LGPL license within the CRAN project

GP-9s Are Ubiquitous Proteins Unlikely Involved in Olfactory Mediation of Social Organization in the Red Imported Fire Ant, Solenopsis invicta

The red imported fire ant (RIFA), Solenopsis invicta, is an invasive species, accidentally introduced in the United States that can cause painful (sometimes life-threatening) stings to human, pets, and livestock. Their colonies have two social forms: monogyne and polygyne that have a single and multiple functional queens, respectively. A major gene (Gp-9), identified as a putative pheromone-binding protein on the basis of a modest amino acid sequence identity, has been suggested to influence the expression of colony social organization. Monogyne queens are reported to possess only the GP-9B alleles, whereas polygyne queens possess both GP-9B and GP-9b. Thus, both social forms are reported to express GP-9B, with GP-9b being a marker expressed in polygynes but it is absent in monogynes. Here, we report two types of polygyne colonies, one that does not express GP-9b (monogyne-like) and the other expressing both proteins, GP-9B and GP-9b. Given their expression pattern, GP-9s are hemolymph proteins, which are more likely to be involved in the transport of lipids and small ligands within the homocoel. GP-9B existed in two forms, one of them is phosphorylated. The helical-rich content of the protein resembles the secondary structures of a beetle hemolymph protein and moth pheromone-binding proteins. An olfactory role is unlikely given the lack of specific expression in the sensillar lymph. In marked contrast to GP-9s, a chemosensory protein, SinvCSP, is demonstrated to be specifically expressed in the antennae. Within the antennae, expression of SinvCSP is restricted to the last two segments, which are known to house olfactory sensilla

Wolbachia and DNA barcoding insects: patterns, potential and problems

Wolbachia is a genus of bacterial endosymbionts that impacts the breeding systems of their hosts. Wolbachia can confuse the patterns of mitochondrial variation, including DNA barcodes, because it influences the pathways through which mitochondria are inherited. We examined the extent to which these endosymbionts are detected in routine DNA barcoding, assessed their impact upon the insect sequence divergence and identification accuracy, and considered the variation present in Wolbachia COI. Using both standard PCR assays (Wolbachia surface coding protein – wsp), and bacterial COI fragments we found evidence of Wolbachia in insect total genomic extracts created for DNA barcoding library construction. When >2 million insect COI trace files were examined on the Barcode of Life Datasystem (BOLD) Wolbachia COI was present in 0.16% of the cases. It is possible to generate Wolbachia COI using standard insect primers; however, that amplicon was never confused with the COI of the host. Wolbachia alleles recovered were predominantly Supergroup A and were broadly distributed geographically and phylogenetically. We conclude that the presence of the Wolbachia DNA in total genomic extracts made from insects is unlikely to compromise the accuracy of the DNA barcode library; in fact, the ability to query this DNA library (the database and the extracts) for endosymbionts is one of the ancillary benefits of such a large scale endeavor – for which we provide several examples. It is our conclusion that regular assays for Wolbachia presence and type can, and should, be adopted by large scale insect barcoding initiatives. While COI is one of the five multi-locus sequence typing (MLST) genes used for categorizing Wolbachia, there is limited overlap with the eukaryotic DNA barcode region

Loss of Wolbachia infection during colonisation in the invasive Argentine ant Linepithema humile

WOLBACHIA are maternally inherited bacteria, which are very common in arthropods and nematodes. Wolbachia infection may affect host reproduction through feminisation, parthenogenesis, male-killing, cytoplasmic incompatibility and increased fecundity. Previous studies showing discrepancies between the phylogenies of Wolbachia and its arthropod hosts indicate that infection is frequently lost, but the causes of symbiont extinction have so far remained elusive. Here, we report data showing that colonisation of new habitats is a possible mechanism leading to the loss of infection. The presence and prevalence of Wolbachia were studied in three native and eight introduced populations of the Argentine ant Linepithema humile. The screening shows that the symbiont is common in the three native L. humile populations analysed. In contrast, Wolbachia was detected in only one of the introduced populations. The loss of infection associated with colonisation of new habitats may result from drift (founder effect) or altered selection pressures in the new habitat. Furthermore, a molecular phylogeny based on sequences of the Wolbachia wsp gene indicates that L. humile has been infected by a single strain. Horizontal transmission of the symbiont may be important in ants as suggested by the sequence similarity of strains in the three genera Linepithema, Acromyrmex, and Solenopsis native from South and Central America

Estimation of Ligament Loading and Anterior Tibial Translation in Healthy and ACL-Deficient Knees During Gait and the Influence of Increasing Tibial Slope Using EMG-Driven Approach

The purpose of this study was to develop a biomechanical model to estimate anterior tibial translation (ATT), anterior shear forces, and ligament loading in the healthy and anterior cruciate ligament (ACL)-deficient knee joint during gait. This model used electromyography (EMG), joint position, and force plate data as inputs to calculate ligament loading during stance phase. First, an EMG-driven model was used to calculate forces for the major muscles crossing the knee joint. The calculated muscle forces were used as inputs to a knee model that incorporated a knee–ligament model in order to solve for ATT and ligament forces. The model took advantage of using EMGs as inputs, and could account for the abnormal muscle activation patterns of ACL-deficient gait. We validated our model by comparing the calculated results with previous in vitro, in vivo, and numerical studies of healthy and ACL-deficient knees, and this gave us confidence on the accuracy of our model calculations. Our model predicted that ATT increased throughout stance phase for the ACL-deficient knee compared with the healthy knee. The medial collateral ligament functioned as the main passive restraint to anterior shear force in the ACL-deficient knee. Although strong co-contraction of knee flexors was found to help restrain ATT in the ACL-deficient knee, it did not counteract the effect of ACL rupture. Posterior inclination angle of the tibial plateau was found to be a crucial parameter in determining knee mechanics, and increasing the tibial slope inclination in our model would increase the resulting ATT and ligament forces in both healthy and ACL-deficient knees

Molecular Variation at a Candidate Gene Implicated in the Regulation of Fire Ant Social Behavior

The fire ant Solenopsis invicta and its close relatives display an important social polymorphism involving differences in colony queen number. Colonies are headed by either a single reproductive queen (monogyne form) or multiple queens (polygyne form). This variation in social organization is associated with variation at the gene Gp-9, with monogyne colonies harboring only B-like allelic variants and polygyne colonies always containing b-like variants as well. We describe naturally occurring variation at Gp-9 in fire ants based on 185 full-length sequences, 136 of which were obtained from S. invicta collected over much of its native range. While there is little overall differentiation between most of the numerous alleles observed, a surprising amount is found in the coding regions of the gene, with such substitutions usually causing amino acid replacements. This elevated coding-region variation may result from a lack of negative selection acting to constrain amino acid replacements over much of the protein, different mutation rates or biases in coding and non-coding sequences, negative selection acting with greater strength on non-coding than coding regions, and/or positive selection acting on the protein. Formal selection analyses provide evidence that the latter force played an important role in the basal b-like lineages coincident with the emergence of polygyny. While our data set reveals considerable paraphyly and polyphyly of S. invicta sequences with respect to those of other fire ant species, the b-like alleles of the socially polymorphic species are monophyletic. An expanded analysis of colonies containing alleles of this clade confirmed the invariant link between their presence and expression of polygyny. Finally, our discovery of several unique alleles bearing various combinations of b-like and B-like codons allows us to conclude that no single b-like residue is completely predictive of polygyne behavior and, thus, potentially causally involved in its expression. Rather, all three typical b-like residues appear to be necessary

Current Status of a Model System: The Gene Gp-9 and Its Association with Social Organization in Fire Ants

The Gp-9 gene in fire ants represents an important model system for studying the evolution of social organization in insects as well as a rich source of information relevant to other major evolutionary topics. An important feature of this system is that polymorphism in social organization is completely associated with allelic variation at Gp-9, such that single-queen colonies (monogyne form) include only inhabitants bearing B-like alleles while multiple-queen colonies (polygyne form) additionally include inhabitants bearing b-like alleles. A recent study of this system by Leal and Ishida (2008) made two major claims, the validity and significance of which we examine here. After reviewing existing literature, analyzing the methods and results of Leal and Ishida (2008), and generating new data from one of their study sites, we conclude that their claim that polygyny can occur in Solenopsis invicta in the U.S.A. in the absence of expression of the b-like allele Gp-9b is unfounded. Moreover, we argue that available information on insect OBPs (the family of proteins to which GP-9 belongs), on the evolutionary/population genetics of Gp-9, and on pheromonal/behavioral control of fire ant colony queen number fails to support their view that GP-9 plays no role in the chemosensory-mediated communication that underpins regulation of social organization. Our analyses lead us to conclude that there are no new reasons to question the existing consensus view of the Gp-9 system outlined in Gotzek and Ross (2007)

Application of Equilibrium Models of Solution Hybridization to Microarray Design and Analysis

Background: The probe percent bound value, calculated using multi-state equilibrium models of solution hybridization, is shown to be useful in understanding the hybridization behavior of microarray probes having 50 nucleotides, with and without mismatches. These longer oligonucleotides are in widespread use on microarrays, but there are few controlled studies of their interactions with mismatched targets compared to 25-mer based platforms. Principal Findings: 50-mer oligonucleotides with centrally placed single, double and triple mismatches were spotted on an array. Over a range of target concentrations it was possible to discriminate binding to perfect matches and mismatches, and the type of mismatch could be predicted accurately in the concentration midrange (100 pM to 200 pM) using solution hybridization modeling methods. These results have implications for microarray design, optimization and analysis methods. Conclusions: Our results highlight the importance of incorporating biophysical factors in both the design and the analysis of microarrays. Use of the probe ‘‘percent bound’ ’ value predicted by equilibrium models of hybridization is confirmed to be important for predicting and interpreting the behavior of long oligonucleotide arrays, as has been shown for shor